Establishing spin-network topologies by repeated projective measurements 
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It has been recently shown that in quantum systems, the complex time evolution typical of many- 
bodied coupled networks can be transformed into a simple, relaxation-like dynamics, by relying on 
periodic dephasings of the off-diagonal density matrix elements. This represents a case of "quantum 
Zeno effects", where unitary evolutions are inhibited by projective measurements. We present here 
a novel exploitation of these effects, by showing that a relaxation-like behaviour is endowed to the 
polarization transfers occurring within a iV-spin coupled network. Experimental implementations 
and coupling constant determinations for complex spin-coupling topologies, are thus demonstrated 
within the field of liquid-state nuclearmagnetic resonance (NMR). 
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Introduction. — A counterintuitive finding in quantum 
theory concerns the temporal evolution of systems sub- 
ject to repeated observations. In the limit of very fre- 
quent measurements, this evolution can be slowed down 
almost indefinitely, in what has been nicknamed the 
quantum Zeno effect [ij. An opposite form can accel- 
erate certain aspects of unitary evolution of a quantum 
system, by relying on other forms of projective measure- 
ments |2j . A particularly suitable field to explore and ex- 
ploit these Zeno and anti-Zeno quantum regimes, arises in 
NMR. This connection was first noted within the frame- 
work of a 'toy' nutation experiment Q. Very recently, 
we have exploited repeated projective measurements in a 
system composed of spins (qubits) entangled with an ef- 
fective lattice [ij], to steer polarization transfers between 
the qubits under mis-matched resonant transfer condi- 
tions into either a quantum Zeno or an anti-Zeno regime. 
In a series of related studies, Tycko used pseudo-random 
timings for switching the coherent character of dipole- 
dipole dynamics to an incoherent limit 

The present study explores the uses of projections 
within the context of Total Correlation Spectroscopy 
(TOCSY), a widespread tool in the arsenal of liquid- 
state NMR Q that serves to establish spin-coupling net- 
work topologies. The TOCSY experiment focuses on a 
homonuclear network of exchange-coupled spins, whose 
inequivalent chemical identities (as defined by the action 
of a chemical shift term) are made equivalent owing to 
the application of a train of suitable, continuous it pulses. 
Under such conditions the effective spin Hamiltonian be- 
comes 
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leading to complex magnetization transfers modes among 
all A^-coupled spin-| partners in a network. As shown in 
Figs. [Tp, QJ on further explained "projections" can be 
used to guide spin evolution to a desired target state in 
an almost monotonic fashion. This in turn can facilitate 
the extraction of information of the Hj coupling struc- 



ture when dealing with complex spin-spin topologies. In 
general, the principles that are here introduced based on 
repeatedly erasing off-diagonal coherences amount to an 
effective switch of an A^-qubit quantum evolution from a 
Hilbert space of dimension 2^, to a simpler, incoherent 
form occurring in an A^-dimensional space. 

Principles and Methods. — In the following we assume 
an N spin-system, subject to the Hamiltonian in Eq. [1] 
Even when containing a small number of coupling con- 
stants {Jij} 1<:i j <N such Hj can generate a complex time 
evolution, owing to its many-body nature Jf|. In fact, 
the multi-spin response originated by the Hamiltonian 
Hj can be exploited for transferring coherences among 
all homonuclear spins within a contiguous structure, and 
thereby help to establish spin-spin connectivities. Fig- 
ure [T^, illustrates the complex kind of spin transfer pat- 
terns that can be promoted by Hj. This complexity can 
eventually become a drawback, as it may demand a sys- 
tematic scanning of evolution times to ensure that no 
spin-spin correlations are fortuitously missed. It may also 
complicate the quantitative determination of the pairwise 
spin-spin coupling parameters. As illustrated by calcu- 
lations (Figs. [lb, [TJ;), projective measurements could 
remedy these disadvantages, and transform the complex 
evolution imposed by H j into a quasi-monotonic transfer. 
The central parameter defining the nature of this transi- 
tion to a Zeno-like relaxation dynamics is the time delay 
r between projective measurements. As this isotropic J- 
evolution period is reduced, the transfer patterns among 
the spins change their originally complex oscillations to 
a smooth and almost monotonic polarization sharing. As 
r becomes very short, this turns into a "freezing" (apart 
from relaxation effects) of all evolutions. It is worthwhile 
noting that for optimally chosen r values nearly equal dis- 
tributions of polarizations - as usually desired in TOCSY 
- can be rapidly achieved. In the absence of instrumental 
losses these amplitudes converge to the ratio of initially 
polarized spins and the total number of sp ins determined 
by the quasi-Boltzmann equilibrium [10]; in Fig. [T]this 
is §=0.4. 
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FIG. 1: (Color online) Time evolutions of single quantum spin-operators h x in Pyridine, (a) Coherent oscillatory transfer 
patterns (simulated — > solid/dashed lines, experimental data — > spin I2, blue triangles, spin I2, green circles, spin ^3, red 
squares) driven by Hj between the spin sites in Pyridine (J-couplings taken from 0]). Implementations involved a selective 
preparation of I^x (a,b) and Ii x (c), followed by the application of a train of n pulses to suppress inter-site chemical shift 
differences. (b,c) Diagrams illustrating the switch of the dynamics shown in (a) to quasi-monotonic polarization transfers, 
as a result of introducing repeated projective measurements. These calculated curves involved instantaneous erasements of 
the off-diagonal density matrix terms at intervals r=31.7ms (optimal transfer), r=12.7ms (intermed.), r=0.63ms (Zeno) for 
panel (b) and r=63.5ms (optimal transfer), r=26.6ms (intermed.), r=3.17ms (Zeno) for panel (c). Points in panel (a) (and 
all remaining data in this work) were acquired on a 600 MHz NMR spectrometer using 50 mM Pyridine in CDCI3. The 
time evolutions of the single quantum spin-operators U x were monitored for each of the three chemically inequivalent sites 
in Pyridine in a point-by-point fashion. The transfer pattern under a free evolution driven by Hj was recreated using pulse 
sequences capable of efficiently suppressing chemical shift differences, including DIPSI-3 §| and MLEV-8 [f| 



Figure [5^ sketches how isotropic mixing sequences can 
be modified to incorporate the projective measurements 
leading to the behaviours in Figs. [Ud, QJ. In lieu of a 
full 2D acquisition where every spin is excited and la- 
beled to follow its dynamics, sequences started with a 
selective pulse followed by a gradient which either ex- 
cite or deplete the signal contributed by an individual 
spin/chemical site. This continues with a looped evolu- 
tion incorporating an isotropic Hj action of length r, a 
storage of the resulting transverse Ii X coherences along 
the z axis of the Bloch sphere, a gradient-based purging 
segment (t p ), and a recall of the stored polarizations for 
further evolutions. This "mixing" section of the sequence 
was looped M times. Central in this sequence is the purg- 
ing segment charged with the erasing of all off-diagonal 
density matrix elements. This consists of a chirped it 
pulse acting in combination with agradient G c hirp for de- 
phasing zero-quantum elements followed by a spoil 
gradient G s . The amplitudes of these gradients were var- 
ied in a random fashion throughout each loop in order to 
avoid fortuitous echoes of the coherences [J]. It is illus- 
trative to examine how the gradient-based elements in 
Fig. [2k help to erase the off-diagonal coherences created 
by the action of Hj. It can be shown that in addition 
to single-spin coherences, this isotropic Hamiltonian will 
lead predominantly to two-spin elements involving zero 
quantum (Ii±Ij T ), single quantum (I iz Ij + ) and double 
quantum (Ii±Ij±) spin-operators @. Even after a ^ stor- 
age pulse the higher-quantum elements in this series will 



possess an off-diagonal character, enabling their cancella- 
tion over the full sample volume when acted upon by the 
G s gradient (Fig. f2b). Such dephasing will not happen 
for the zero quantum terms, yet these will disappear by 
the action of the frequency-swept pulse applied in com- 
bination with a random field gradient (Fig. |2t). 

With this sequence at hand, the experimental imple- 
mentation of the "projected" Hj evolution was assayed 
on Pyridine's protons for a number of conditions. A se- 
ries of experimental data are displayed in Figs. [3^ and 
|3Jd, showing the redistribution of polarization following 
the selective excitation of sites Ix/Iy and h/hh respec- 
tively. The behaviour of each site's magnetization is fun- 
damentally different from that shown in Fig. QJ, under 
the action of Hj, but agrees well with the expectations 
presented in Figs. [TJd, [TJ;. The main discrepancies be- 
tween the experimental and the simulated data arise from 
the fact that the asymptotic polarization distribution is 
reached at a lower signal amplitude (~0.3) than ideally 
expected (0.4); rf pulse inaccuracies and/or relaxation- 
derived losses are probably responsible for this. A related 
switching to a relaxation-like evolution can be observed 
if the selective excitation of a particular site, is replaced 
by its selective demagnetization. The corresponding ex- 
perimental data (Fig. [3};), illustrates the possibility of 
rcpolarizing a depleted magnetization via neighbouring 
J-coupled spins. For the Pyridine system assayed, the 
theoretical signal amplitudes that should be achieved in 
such equalization experiments are given by | =0.6 - very 
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FIG. 2: (Color online) (a) Polarization transfer sequence imitating projective measurements interspersed with segments of 
isotropic Hj evolution. Commonly used parameters were M=20-40, £ se ;=5-10ms, r=l-60ms, t p =25-30ms, G s =20 G/cm, 
G c hirp~0-5-l G/cm and a sweep range of the chirped n pulses of Ao;=20kHz. (b,c) Transfer characteristics expected during 
the course of this sequence for one-spin single quantum (b) and two-spin zero and double quantum spin-operators (c) . Notice 
the predicted suppressions of off-diagonal elements of the density matrix over the full sample volume, and the steering of the 
(Iim) transfer patterns towards a monotonic evolution. 



close to the experimentally obtained values. 

It follows that projective processes can transform a 
complex dynamics like that driven by Hj, into a quasi- 
monotonic behaviour. It is worth exploring to what ex- 
tent can the coupling constants effecting these transfers, 
be extracted from such relaxation-like dynamics. This is 
facilitated by the short-r regime leading to the kind of 
curves illustrated in Fig. |31 where the original 2 N dimen- 
sionality defining the coupled spins' wavefunctions is pro- 
jected into an ./V-dimensional space. In such a space only 
single spin polarizations do matter. The corresponding 
spin dynamics can then be represented by a polarization 
vector P(M r) 



P(M r) = [p(r)] M P(0). 



(2) 



This leads, after M projections, to a two-spin polariza- 
tion matrix 
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The simple quadratic behaviour evidenced by the off- 
diagonal build-up in this matrix, can be expanded to 
the many-body Hamiltonian of Eq. [TJ Assuming then 
that the condition rJij<Cl still holds for all pairwise cou- 
plings, the on- and off-diagonal elements of [p(r)) M can 
be approximated by 
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Here M depicts the number and r the timing of the pro- 
jective measurements. The vector P(0) is the initial po- 
larization of each spin i, and the polarization transfer ma- 
trix [p{t)] m represents an Nx N, two-dimensional spec- 
trum whose elements carry the intensities of the self- and 
of the cross-correlations observed at times M ■ r. These 
intensities can be simply computed from the knowledge 
that under the action of Hj=JIi ■ Ij, the two-spin evo- 
lution is of the form 0] 
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where O represents higher-order relayed transfer terms. 
Considering short durations r and the fact that projec- 
tive measurements will only preserve diagonal terms, jus- 
tify the approximations 
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This matrix defines the signal build-up of initially un- 
polarized spins i. Therefore, the absolute values of the 
involved couplings can be extracted either from the initial 
build-up crosspeaks observed in a two-dimensional spec- 
tral distribution, from the ratio between peak intensities 
obtained in ID build-up transfer curves, or from these 
crosspeak values normalized by the self-peak intensity. 
For Pyridine - an effective five spin-system - this would 
mean that J13 and an average of J12 and Jy% could be ex- 
tracted from the short-r slopes of the data obtained after 
a selective excitation of I\ (Fig. [3t), etc. Alternatively, 
the coupling constants involved can be established from 
the decaying signal amplitude of the initially polarized 
spin species. For the model system hereby examined the 
individual coupling constants were quantified from the 
ratios of the signals following projective measurements. 
The resulting experimental values are summarized in Ta- 
ble HJ and show a good agreement with values found in 
literature, as well as when compared against the predic- 
tions of numerical simulations. 
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FIG. 3: (Color online) Experimental single spin polarizations, obtained from the sequence in figure [2^, for Pyridine. (a,b) 
Transfer characteristics resulting from a "projected" Hj evolution, following selective preparations of the spins ii (a) and I2 (b) 
for r=63.5ms (filled symbols) and r=12.5ms (empty symbols), (c) Time evolution for an analogous repolarization experiment 
performed after a selective demagnetization of spin Zi for r=50.45ms (filled symbols) and r=I2.61ms (empty symbols). 
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TABLE I: Experimental J- values extracted from short-r peak 
ratios compared against simulated numerical expectations, 
and literature values[7j. 



Conclusion and Outlook. — This study outlined a 
strategy that uses well-known elements from NMR spec- 
troscopy's toolbox, like field gradients and rf pulses, 
to imitate projective measurements. This allowed us 
to use NMR as a "quantum simulator" for describ- 
ing novel Zeno-likc effects in exchange-coupled spin net- 
works. These concepts also lead to a new experimental 
approach with potentially useful applications, whereby 
complex and a priori unknown polarization transfer pat- 
terns are morphed into smooth, monotonic relaxation- 
like functions. From these it is straightforward to de- 
termine the underlying spin-spin coupling constants, to 
establish unambiguous correlations among spin networks, 
or to redistribute pools of unused polarization among the 
spins. Future work will show extensions of this strat- 
egy to additional NMR experiments, including general- 
izations to other types of scalar, dipolar and quadrupo- 
lar multi-spin dynamics. We have also found this ap- 
proach useful to shorten recycle delays - particularly in 
heteronuclear polarization transfer experiments [lS^ . It 
is conceivable that the use of dynamical decoupling via 
frequent 7r-pulses can bring about analogous results |13| . 
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